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The single-file problem of A'' particles in one spatial dimension is analyzed, when each 

particle has a randomly distributed diffusion constant D sampled in a density p{D). The averaged 

one-particle distributions of the edge particles and the asymptotic (A'^ ^ 1) behaviours of their 

transport coefficients (anomalous velocity and diffusion constant) are strongly dependent on the 

^^ , D-distribution law, broad or narrow. When p is exponential, it is shown that the average one- 

^^ ■ particle front for the edge particles does not shrink when A'^ becomes very large, as contrasted to the 

\^^ ' pure (non-disordered) case. In addition, when p is a broad law, the same occurs for the averaged 

front, which can even have infinite mean and variance. On the other hand, it is shown that the 

central particle, dynamically trapped by all others as it is, follows a narrow distribution, which is a 

Gaussian (with a diffusion constant scaling as A'^"^) when the fractional moment < D~^" > exists 

flj ■ and is finite; otherwise (p(-D) oc D"~^, a < ^), this density is, far from the origin, a stretched 

Cj ' exponential with an exponent in the range ]0, 2]; then the effective diffusion constant scales as N~" , 

^ ' with /3 = 1/(2q). 
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TJ . I. INTRODUCTION 

a. 

The single-file diffusion problem is encountered in various fields (one-dimensional hopping conductivity |]| , ion 
transport in biological membranes |g,^, channelling in zeolithes |j]). Generally speaking, this is modelled by a set of 
N diffusing particles on the line with hard-core repulsion; due to such an interaction, any initial ordering is preserved 
\^J ■ in the course of time and the particles can be, once for all, labelled 1, 2, . . . , N from left to right. As in ref Q, I 
^vj vifill here consider the case where, at the initial time, these particles form a compact cluster centered at the origin; 

^^ ' the solution for an arbitrary initial condition has been given in ref [pi, using the reflection principle. In addition, as 
,—1 I contrasted to the standard model in which all the particles have the same diffusion constant D, it is here assumed 
CN| ■ that the particle i has a random Di, chosen independently of all the others in a given distribution law p{D). The 
^D \ latter can always be written as follows: 



where Dq denotes a specific value of the diffusion constant (e. g. the average value) and r(^) a positive function 
normalized to unity: 



!,.._ ^....±- .:: 

K^ , relation, it can result from a random radius; more directly, this also happens for particles having random masses. 
Vh The main purpose of this paper is to analyze how the choice of the distribution r modifies the large- A^ depen- 

5r ! dence of the transport coefficients for one particle of the cluster and, more generally, to find out the density probability 
of its coordinate in the presence of random diffusion constants. More specifically, focus will be given on edge particles 
and on the one which (for A^ odd) is at the middle of the cluster; for the "pure" case (all the particles have the same 
D), asymptotic laws (A^ ^ 1) have been given in ref. H. Among other results, it was shown that increasing A^ yields 
a narrowing of the one particle-probability distributions, as a result of the "pressure" exerted by other particles on a 



'e-mail: aslangul@gps.jussieu.fr 



given one. This will be shown not to be true in all cases when the D's are randomly sampled. As a rule, if the middle 
particle is most often insensitive to disorder - provided that the average < D~^'^ > exists and is finite - it will be 
shown that the distribution for the side particles strongly depends on the latter, going from narrow to broad laws. 



II. BASIC RELATIONS 



The diffusion equation here writes: 
d 



N 



dt 



p{xi,X2,.- ■,XN]t) = ^ Dn -—yP{xi,X2,.- ■,XN;t) 



dxl 



(2.1) 



As a consequence of the hard-core repulsion and of the above-mentionned initial condition, the solution has the 
expression: 



p{xi,X2,-. .,XN 



N -^l/(4D„t) ^-1 
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(2.2) 



In (2.2), Y is the Heaviside unit step function (Y{x) = 1 if a; > 0, otherwise), D„ denotes the diffusion constant of 
the n^ particle and Cn is the normalization constant. Since the D's are assumed to be statistically independent, the 
various averages can be readily performed and the average Af-particle front turns out to be: 



N 



N-1 



<p{xi,X2,...,XN;t)>= NlYl <G{Xn,t)> Y[Y{Xn+l 



(2.3) 



r! = l 



where < G{xn, t) > is the average of the gaussian distribution with the density p(D): 
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/ is an even positive function normalized to ^ on the interval [0, +00]. 

From (|2.3|), all the one-particle averaged fronts can be formally obtained by integrating over all a;'s but one: 
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<p«(x„;t)>= ( n 

I TTi— 1, m^n 



dXm \ <p{xi,X2,-.-,XN\t)> 



One then finds: 



<pW(a^;i)> = 
where I{u) is defined by: 

I{u) = 2 
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(2.6) 



(2.7) 



(2.8) 



where $ denotes the probabihty integral 0. In (p?7|), the two factors [1 ± /(u)]/2 represent the steric effects on the 
n*'^ particle due to all other ones. I{u) is a non-decreasing odd function such that 7(±oo) — ±1 (for the pure case, 
one simply has I{u) = ^{u)). For the edge particles (n = 1 or n = A^), < p„ ix;t) > can be given a form occurring 
also in the theory of extreme events B , namely: 
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<pY'j^{x;t)>= ± 



{4Doty/^ du 
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(2.9) 



The — sign (resp +) refers to the left particle (n = 1) (resp. right particle {n = N)). In any case, the following 
relation holds true: 
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For n arbitrary, (|2.9|) generalizes into: 



<pl^>{x;t)>: 



<pY'{x;t)> = <p'j^>{-x;t)> 



_d_ fl + I{u) 
du 
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For the particle which is at the center of the cluster {N odd), one has: 
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< P{N+i)/2i^: i)>- (4^gi)i/2 {[(Ar_i)/2]!}2 



1-P{u) 
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or, equivalently, from ( ^.ll| ): 
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(2.10) 



(2.11) 
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In the following sections, the above equations will be used with different choices for the distribution r(^). Exact 
asymptotic results for large N can be obtained due to the fact that, for A^ ^ 1, the steric factors {^[1 ± /(m)]}" 
display a rather sharp variation. 

Generally speaking, the one-particle transport coefficients can be found from the behaviour at small k of the 
characteristic function Iin{K): 



/+00 
dx e"'^ < p(,i) (x; t)> {K = k^/AD^t) 

-OC 



(2.14) 



The transport properties of the n"* particle subjected to the random field of all the other are essentially described by 
the two first cumulants: 
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Xn>{t) = Vi/2,niN)Vi 



A: 



<xl> (t) - [< xn > (t)]2 = 2D„{N) t 



(2.15) 



The above time-dependences come from the fact that the only available length-scale is (Doi)^^^, so that < x^ > oc t^/"^ 
at all times and for any A. Thus, at any time, the drift is always anomalous (no velocity) and the mean square 
displacement always has a purely diffusive motion. The following relations are trivially veridied: 



<xi> {t) = - <X]s > (t) 



Axl{t) = Ax^(t) 



(2.16) 



In addition to the one-particle transport coefficients, statistical correlations between the particles can be ana- 
lyzed. For definiteness, focus will be given on the correlations between the two edge particles, altogether contained 
in the two-body density: 
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Using (2.3) - (2.5), one finds 



<pfU=^i,XN;t)>= -±-N{N-l) 
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(2.18) 



III. EXPONENTIAL DISTRIBUTION 



As a first simple example, let us choose: 
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Here, Dg is the expectation value of D. From (2.5) and (|2.8|), one readily obtains 
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I{u) = sgnu(l-e"2|"l) 



The reduced densities for the edge particles density can now be explicitely calculated from (2.9) with: 
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(3.1) 



(3.2) 



(3.3) 



and P+{u) — P_(— li). From (3^), the generating functions IIi^n{K) ( 2.14 ) for the edge particles are found as the 
following: 
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Now, by expanding Hi. jv in the vicinity oi K = 0, and by coming back to the x-variable, one obtains: 

<XN>it) = - <xi> (t) = [ln(A^/2) + C + 0(2-^)] ^^t = yi/2,edge(^) V^ 
where C is the Euler's constant (C = 0.577. . .). The mean square displacements are: 



Ax% 



Aa;? = [y + 0(2-'')] Dot = 2Dedge(iV)t 



(3.4) 



(3.5) 



(3.6) 



As compared to the pure case, the large- A^ behaviour of the transport coefficients is frankly different. When all the 
particles have the same diffusion constant, one has 0: 



nS.e(^) « (1-^)^ 



D^^:\n) oc (InA^)- 



(3.7) 



As contrasted, for exponentially distributed random D, the anomalous drift coefficient increases like IniV, 
whereas the effective diffusion constant D^^„q{N) tends towards a finite (non- vanishing) value. Also note that the 
height of the maximum also saturates and does not go to zero for N infinite. As a whole, the effect of an exponential 
disorder is rather subbtle: most particles of the cluster have a rather small D, so that the pressure they exert on the 
edge ones is smaller than in the pure case; on the other hand, large D forces the packet to more with a "velocity" 
which increases more rapidly with N than in the pure case. The interplay of these facts produces a distribution 
having a finite height and width at infinite N. This is illustrated on Fig. ||, which clearly shows the width and height 
saturation when TV becomes very large, as quantitatively described in eq. (^.6|). It is readily seen that the value of 



^ Pi N ^ at its maximum is close to (2/e) oc A^". 

For A^ ^ 1, the abscissa of the edge particles are approximately distributed according to: 



<pW(x;t)>~ Y{u) 
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<pY>{x;t)>^<p'j^>{~x;t)> 
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(3.8) 



These functions are exponentially small (~ e ^) near | m | = and display a plain exponential decay (^ e ^'"1) for 

I M |» ln[(iV/2)^/^]. < p}^' > is maximum for u — Umax = ^ ln(A^/2); this shows that the corresponding Xmax 
coincides with ± < xjy > (see (3^)). The approximate expressions (3^) are hardly distinguishable from their exact 
counterparts as soon as N is greater than a rather small number (~ 50) and can be rewritten in terms of the proper 
shifted (but here not rescaled) variable X: 



<p«(x;i)>^ Y{x) 



{Doty/^ 
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(3.9) 



The statistical correlations between the two edge particles are most simply measured by the correlator Cin- 

CiN — < xi xn > — < xi > < xn > (3.10) 



which, due to scaling in space, varies hnearly in time. Cin can be found by using ( 2.18 ) and (3^); for A^ ^ 1, a httle 
algebra yields: 



CiN - 4Dot [ln^(Ar/2) + 0(lniV)] 



(3.11) 



By (3.(:), this implies that the normalized ratio Ciat/Ax^ ^ - which is constant in time - has a logarithmic increase 
with N: 
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(3.12) 
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Although the numerator and the denominator of the ratio have separately different behaviours as compared to the 
pure case, the ratio still has a ~ (In A^)^ increase with N, exactly as in this latt er case 0. 
For the particle located at the middle of the cluster, the expression (2.13) writes: 
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From (3.13), the following large- iV expression of the mean square dispersion is obtained: 
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where B denotes the Euler beta-function m. Expansion for large- Af yields: 



Ax?^,^1^/^(a;;i) ~ tt-Do^ 
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and provides the asymptotic dependence of the diffusion constant: 

^middle (^^) - ^^0 
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This behaviour is the same as for the pure case S. Thus, for the exponential distribution, the transport coefficient 
of the central particle is essentially unaltered by disorder - as shown below this is also true for any r density having 
< D^^/'^ > finite. The asymptotic density of the middle particle is obtained from (p. 13) as the following: 
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(3.17) 



This function is normalized to unity up to A^^^ terms. Apart from the cusp at u = and outside the region | w | < l/A'^, 
this is essentially a gaussian distribution, exactly as in the pure case. For large A^, the cusp can even be ignored and 
the distribution looks like a Gaussian everywhere. 

In fact, it can now be stated the conditions for the stochastic dynamics of the central particle to be unaffected 
by disorder; such a particle is stymied in some way due to other particles and it can be suspected that strong steric 
effects most often dominate, except perhaps when the distribution p{D) gives a hi gh p robability to find very small D. 
The analysis goes as follows. For A^ ^ 1, it is readily seen that the distribution (2. IS) ass umes non-negligible values 
only for \u |<C 1. This in turn allows to replace I{u) by its small-M expansion; from (2.8) and (|]j), /(O) exists and 
is finite provided that p{D) is bounded by D~^ (/3 < i) for D ^ 0. This means that when p diverges at small D, 
this divergence is not too severe, which implies that finding a small D has not such a great probability. With this 
assumption, one has: 
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I{u) ^ 2uf{0) /(O) ^- J^ da-'^'riO ^ ^ < (^^j > . (3.18) 

Using Stirling formula, and expanding e^'^^^^^/^l in[i-4« / (o)]^ ^^^ eventually gets the asymptotic approximation: 
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<pg,,/,(.;i)>. (^^j /(0)e--/^(«)"\ (3.19) 

Note that the approximation automatically generates a properly normalized density. So, for any p{D) such that 
lim£)^o[-D^'^p(-D)] == 0, the middle particle is essentially distributed according a normal law, the width of which 
decreasing as N^^''^, up to unessential numerical factors. The expression ( |3.19 ) can be written as a universal law in 
term of the proper scaled variable X: 



expressing as a whole the irrelevance of details about the D distributions for the central particle when p{D) is bounded 
near I? = as stated above, which entails that the fractional moment < D~^/^ > exists and is finite. The opposite 
case is treated in the next section which deals with the Gamma-distribution. 

IV. GAMMA DISTRIBUTION 

As another example - which in fact contains the exponential distribution as a particular case - let us choose a 
Gamma distribution: 

PiD) = p(^ D'^-' e^^/^» (a > 0) . (4.1) 



With (|4.1| ), the expectation values of D and D are respectively equal to aDa and a{a + 1)Dq. The exponential 



distribution is recovered by sett ing a — 1, whereas the pure case can be obtained by taking the limit a — > oo. Do — > 0, 



aDo = const. The / function (2.5) is 



f{u) = ^|— \ur^ K i(2|u|) (4.2) 

where K i is the Bessel function of imaginary argument M]. /(O) exists and is finite if a > ^■ 
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Vtt (a- ^)T{a) 2 

The large- A^ approximation of the one-particle densities can be obtained along the same lines as in the previous 
section, by using the approximate expression: 

/(u) ^l^J-\u r-i e-2|"l (I u |» 1) . (4.4) 

T{a) 

Let us consider first the right particle; one has: 

<pW(.;t)>. n")(^rRe--^---%-e- . (4.5) 

When a < 1, i. e. when p diverges at small D, < pjy > is exponentially small near a; = 0. On the contrary, for a > 1, 
< p]y > strictly vanishes as u°'~^, but this behaviour is realized on a very small interval, namely u < N^^'^"'~^\ 
beyond which < p)^ > remains exponentially small. On the other hand, at large u, < p)^ > has essentially an 
exponential decay, so that all the average values < a;™ > exist and are finite. < p]^ > is maximum at Mmax: 

1 N 



Again, Xmax = (4-Do^)^^^ "max is the relevant rescaled variable, since (4.5) can be rewritten as follows: 



<p«(x;t) >^ r(u)^^^i^e-("/"-)°-^-^'"^"-' (V-max)"-^e-2("/"-) _ (47) 



This readily gives < p)^ >maxfx N° which in turns allows to guess that, as for the exponential case, one has 
< a:jv >oc InA^ and < Axf ^ > oc 7V°; it thus turns out that detailed calculations for a ^ 1 are indeed of little 
interest. Up to unessential numerical factors, the large-A'' asymptotic laws all belong to the same class for Gamma- 
distributed diffusion constants, the behaviours being rather easily obtained by considering the simpler exponential 

case a — I; the same obviously holds for statistical correlations between edge particles. 

For the central particle, one has to distinguish the two cases a > ^ and a < ^. Due to (|4.3|), the first one is 



described by the general expression (3.2C). On the contrary, for a < ^j one has: 
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This in turns entails that: 
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As a consequence, one obtains: 
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(4.8) 



(4.9) 



(4.10) 



This shows that when the probability is very high to find quite small diffusion constants, the distribution for 
the central particle diverges at a; = (but is clearly integrable) and is essentially a stretched exponential when | u \ 
is large: 
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For a = i /(m) is proportionnal to the Bessel function Kq and the divergence is logarithmic. In this case, one has: 



< P(Ar+i)/2l2;, t)>- p^^^i/2 V ^3 ° ^'^ I " I 
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which entails: 
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The divergence at a; = com es from the interplay of pressure and of the frequent occurrence of central particle with 
a very small D. From ( 4.10 ), one finds: 
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Summing up, the effective diffusion constant for the central particle and for the Gamma distribution scales 
with N as follows: 
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V. BROAD DISTRIBUTIONS 



The two previously examples discard the interesting case in which the diffusion constants are sampled in a broad 
law, possibly devoid of usual first few moments (mean and variance). Obviously enough, such a case is interesting 
since the average front for the edge particles are expected to be strongly asymmetric and also quite diffuse, due to 
the pressure exerted by inner particles and the possibility of rather high diffusion contants. This fact even opens 
the possibility for the edge particules to be also distributed according a broad law, devoid of mean and variance. 
Generally speaking, I will now consider the case: 
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All the moments < D^ > div erge for fc > /i 
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With the distribution (5.1), the / function (2.4) can be expressed in terms of the incomplete 7 function 
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Again, I am interested in the large- iV limit, in which case the right particle density (2.£) assumes non- negligible values 
for u of the order of or greater than uq defined as: 
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For genuine broad laws {fi small), uq is much greater than 1 even for moderately large values of N. This allows to state 
again that the large-u expansion of I{u) is only relevant, and to substitute everywhere the approximate expression: 
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Injecting this in ( |2.9| ) yields the large-iV approximation, valid for u > 0: 

< pW (:,; ^) > c. iV ^ e- ^ r(,+ i)/(2.v. „..) ^_^,^^,^ 7(^+1,^.2) 
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it being understood that < p)^ > is exponentially small for u < and can be considered as identically vanishing for 

0, namely: 



M < 0; note that the expression (5.5) goes toward zero extremely rapidely when u 

<p«>^e-Cst^/«^- 
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On the other hand, this distribution is indeed a broad law in the wide sense, since (5^) displays for any /i a power-law 
behaviour at (very) large u: 



<p^^\x;t)>~NJ^T{ti+^)u-^^'^+'^ (w»7Vi/(2A')) 
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Eqs. ( p.6[ ) and ( p.7\ ) show that the average front for the right particle is strongly asymmetric, displaying a rather steep 
increase on the left (towar ds th e inner part of the cluster) and a very slow decrease on the other side, towards the free 
part of space. The result ( |5.7| ) entails that the moment < u™ > exists and is finite only if the following inequality is 
satisfied: 

m < 2/i . (5.8) 

Thus, for n < h, the expectation value of x (as well as all higher moments) is infinite. The one-particle distribution 
(plotted in Fig. for a few values of A'^) can nevertheless be characterized by the value Mmax giving its maximum value, 
which turns out to coincide with uq (see ( ^.3[ )), except for numerical factors. It is readily verified that, accordingly, 
the value of < p-j^ ^ > at its maximum is ex A^~ i/(2m)_ As a whole, the maximum of < pj^ {x;t) > moves in time as 
the following: 
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Xjxiax is here the relevant rescaled coordinate; agreement is quite good with the exact results (see Fig. 0). 



n(l) 



Note that for /i ^ 00, (5.7) shows that < p)^ > decreases faster than any arbitrary power of u , in agreement 



with B . Also note that no obvious measure of correlations here exists since for /i < 1 the mean square dispersions 
are infinite. 



VI. SUMMARY AND CONCLUSIONS 

For the single-file diffusion problem with random diffusion constants, the asymptotic laws (large N) giving the 
transport coefficients and the averaged one-particle densities have been derived. 



Generally speaking, the asymptotic distribution of the central particle is most often gaussian (see (3.19)), as 
a result of the fact that steric effects nearly always dominate; when this is the case, the dynamical trapping of this 
particle is insensitive to details describing the distribution p{D) of the diffusion constants. On the other hand, when 
p{D) diverges as D~^^^ or faster at small D {p{D) ex 13"^^, a < ^), the distribution of the central particle is no 

more a no rmal one: at small x, it itself goes to infinity as a;~*^^~^"^ and is a stretched exponential e~'^^ " at large x 
(see ( 4.10| )). The effective diffusion constant scales as A^~i/(2«) for a < ^ and as A''"-^ otherwise, showing that in the 
iV- infinite limit, as in the pure case, the motion indeed becomes subdiffusive in agreement with Harris result [^. 

On the contrary, the stochastic dynamics of the edge particles is strongly dependent upon the nature (narrow 
or broad) of the diffusion constants probability density p{D). When the latter is narrow (i. e. exponential), the 
effective diffusion constant of edge particles tends toward a finite constant when N becomes infinite; the average 
coordinate increases with N as InA^, faster than for the pure case. When the p{D) is broad (i e. behaves as a 
power law £)~(''+i) at large D), the edge particles are themselves distributed according a broad law, devoid of mean 
and variance when p is smaller than i. The distribution of the coordinate, strongly asymmetric (see Fig. 0), can 



nevertheless be characterized by the abscissa of its maximum (5.9); the prefactor of the latter displays a A^i/f^A") 
scaling i. e. rapidly increases with N when p is small. This comes from the fact that there is a high probability to 
find inner particles with a large diffusion constant which, in turns, gives rise to a high "pressure" exerted by the core 
of the cluster on its "surface" . 

In all cases, the averaged one-particle asymptotic density has been found; it can be written under a general 
form displaying the basic ingredients of the problem. For the right particle, one has: 



Ji) 



<p-j^'{x] i) >~ iVexp 



N 



< 1 - $[a;/(4Dt)i/2] > 



< G{x, t) > 



(6.1) 



where < . . . > denotes averaging with the distribution p{D) of the diffusion constants, <& is the probability integral 
and G is the gaussian distribution. For the central particle {N odd), one has: 



< 



„(i) 



P(N+l)/2\^'^ *) -*' 



N 
2^t 



1/2 



< 



15-1/2 



> cxp 



-A^< D-^'^ > 



27rt 



(6.2) 



provided that the fractional moment < D ^/^ > exists and is finite. Otherwise, when p ex D" ^ at small D with 
a < |, < P(N+i)/2 > i^ "^f the form (see ( 1.10 ) for details): 



< pS + 1)/2(^; t)>^C \x r(l-2a) g^p [_ C' I X I 



4al 



(6.3) 



At large x, this gives a stretched exponential with an exponent in the range ]0, 2], whereas the density diverges at 
a; = as a power-law (ex x~^'^°'~^'); for the particular case a = ^, this singularity becomes logarithmic. 
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Figure Captions 



J^)i 



1. Exact averaged front < Pi {x;t) > for the left particle with exponentially distributed diffusion constants; the 
abscissa is the dimensionless variable u = x/{4DQty^^. Each curve is labelled by the number of particles in the 
cluster. 



.(1) 



Averaged front < P(7v+i~i/2('''' ^) ^ '-'^ ^^^ central particle (4.10) for a cluster of 10 particles when the dif- 
fusion constants are distributed according to (4J.) with a = j, the abscissa is the dimensionless variable 

3. Exact averaged front < p}^ {x; t) > for the right particle when the diffusion constants of the cluster arc distributed 
according to (O) with ^i — ^] the abscissa is the dimensionless variable u = a;/(4£'oi)^'^- Each curve is labelled 
by the number of particles in the cluster. 
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